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Abstract 

This  paper  is  concerned  with  the  dynamics  of  an  elastic  bar  that  can  undergo 
reversible  stress-induced  phase  transformations.  We  consider  a  Riemann  problem 
in  which  the  initiaJ  strains  belong  to  a  single  metastable  phase  and  prove  uniqueness 
of  solution  that  satisfy  a  nucleation  criterion  and  a  kinetic  law  at  all  subsonic  and 
sonic  phase  boundaries.  This  paper  generalizes  the  results  of  [3];  the  authors  of  [3] 
considered  a  piecewise  linear  material  for  which  no  wave  fans  exist,  shock  waves 
always  travel  at  the  acoustic  speed  and  shock  waves  are  dissipation-free.  The 
material  model  of  the  present  paper  does  not  suffer  from  these  degeneracies. 

1.  Introduction.  In  the  simplest  one-dimensional  theory  describing  the  longitudi¬ 
nal  motions  of  an  elastic  bar.  one  employs  a  pair  of  conservation  laws  associated  with 
momentum  balance  and  kinematic  compatibility.  When  the  motion  of  the  bar  involves  a 
propagating  strain  discontinuity,  it  is  subject  to  a  pair  of  jump  conditions  associated  with 
these  conservation  laws.  In  addition,  the  second  law  of  thermodynamics  requires  that  the 
dissipation  associated  with  the  moving  discontinuity  be  non-negative,  a  condition  usually 
referred  to  as  the  entropy  inequality. 

The  character  of  the  material  of  the  bar  enters  the  conservation  laws  and  jump  con¬ 
ditions  through  the  stress-strain  relation  a  =  <7(7).  If  the  material  of  the  bar  is  such  that 
stress  is  a  monotonically  increasing  function  of  strain  that  is  strictly  convex  or  strictly 
concave,  then  phase  transformations  cannot  occur,  and  all  propagating  discontinuities 
are  shock  waves.  For  a  bar  made  of  such  a  material,  it  follows  from  a  result  of  Oleinik 
[14]  that  the  Cauchy  problem  for  the  associated  field  equations  and  jump  conditions  has 


1 


at  most  one  piecewise  smooth  solution  that  fulfills  the  entropy  inequality;  see  [13]  for  a 
discussion  of  Oleinik’s  theorem  and  related  results. 

When  the  material  can  undergo  a  reversible,  or  thermoelastic,  phase  transformation, 
the  one-dimensional  elastic  continuum  can  be  characterized  by  a  non-monotonic  relation 
(7  =  <5’(7)  between  stress  and  strain,  or  equivalently  by  a  nonconvex  elaistic  potential 
Typically  one  encounters  a  stress  response  function  <7(7)  in  which  stress  first 
increases  with  increasing  strain,  then  decreases,  and  finally  increases  again;  for  example, 
Ericksen  [10],  The  rising  branches  of  such  a  stress-strain  curve  are  identified  with  different 
phases  of  the  material,  while  the  declining  branch  is  associated  with  an  “unstable  phase.’’ 
For  suitable  value  of  stress,  the  associated  potential  energy  Gi'),a)  =  H'(7)  —  cry  has 
multiple  energy-wells,  each  energy-well  being  associated  with  a  distinct  phase  of  the 
material.  During  a  typical  thermomechanical  process,  the  material  often  moves  from 
one  energy-well  to  another,  or  equivalently,  from  one  branch  of  the  stress-strain  curve  to 
another. 

When  the  stress-strain  curve  is  non-monotonic  and  undergoes  a  change  in  the  sign  of 
its  curvature,  the  Cauchy  problem  need  no  longer  have  a  unique  solution,  even  with  the 
entropy  inequality  in  force;  see  the  remarks  of  Dafermos  [9].  In  order  to  secure  unique¬ 
ness,  many  researchers  have  replaced  the  entropy  inequality  with  various  “admissibility 
conditions”  which  are  to  be  satisfied  by  the  weak  solutions.  For  example,  two  different 
notions  of  maximum  entropy  production  have  been  proposed  by  Dafermos  [8],  and  aug¬ 
menting  the  elastic  theory  with  viscosity  and  capillarity  effects  has  been  proposed  by 
Slemorod  [15],  Truskinovsky  [17].  The  implications  of  these  criteria  for  dynamic  phase 
transitions  have  been  examined  by,  for  example,  Hattori  [11],  James  [12]  and  Shearer 


A  completely  different  approach  has  been  proposed  and  studied  in  [2-5]  based  on  the 
observation  in  [1]  that  the  lack  of  uniqueness  arises  not  only  in  dynamic  motions,  but 
in  quasi-static  motions  as  well.  It  was  suggested  in  [2]  that  in  addition  to  the  usual  - 

constitutive  law  between  stress  and  strain,  further  material  description  in  the  form  of  a 


_ _ n  ; .^.i;  1  i  ity  7o{j8a 

Av  y  01* 

Dlat  I  Cponial 


han 


nucleation  criterion  and  a  kinetic  relation  pertaining  to  the  phase  transitions  are  needed. 
The  importance  of  a  nucleation  criterion  and  a  kinetic  relation  in  the  description  of 
phase  transitions  in  solids  has  long  been  recognized  in  the  materials  science  literature, 
e.g.  Christian  [7]. 

It  was  shown  in  [2]  that  the  inclusion  in  the  continuum  theory  of  the  nucleation  crite¬ 
rion  and  the  kinetic  relation  leads  to  a  determinate  quasi-static  theory  whose  predictions 
are  in  qualitative  accord  with  experiments  on  shape  memory  alloys  that  involve  s!ovvl\- 
propagating  phase  boundaries.  A  similar  result  in  the  dynamical  setting  for  the  Rie- 
mann  problem  for  a  special  piecewise  linear  elastic  material  was  established  in  [3].  It  was 
shown  in  [4,5]  that  the  maximum  entropy  rate  admissibility  criterion  and  the  viscosity- 
capillarity  admissibility  criterion  may  in  fact  be  viewed  as  being  two  particular  examples 
of  kinetic  relations. 

The  study  in  [3]  was  restricted  to  a  piecewise  linear  elaistic  material.  The  nature  of 
this  trilinear  material  model  leads  to  a  considerable  simplification  in  the  analysis.  First, 
this  material  does  not  sustain  wave  fans.  Second,  shock  waves  always  travel  at  the  sound 
speed.  Third,  shock  waves  are  dissipation-free.  For  these  reasons,  it  is  natural  to  question 
whether  the  results  found  in  [3]  were  special  to  the  trilinear  material  and  inquire  whether 
its  conclusions  hold  for  more  general  rising-falling-rising  stress-strain  curves.  Thi^  is  the 
objective  of  the  present  paper.  We  show  for  a  material  whose  rising-falling-rising  stress- 
strain  curve  is  smooth  and  has  a  single  inflection  point,  that  the  nucleatif  ji  criterion  and 
the  kinetic  relation  (applied  to  all  subsonic  and  sonic  phase  boundaries)  serve  to  single 
out  a  unique  solution  to  the  Riemann  problem  with  initial  data  m  a  single,  melastable 
phase.  If  the  kinetic  law  is  applied  only  to  phase  boundaries  that  are  subsonic  and 
not  to  those  that  are  sonic,  we  do  not  have  uniqueness  of  solution  for  all  initial  data. 
Truskinovsky  [18]  takes  the  view  that  the  kinetic  relation  should  not  be  applied  to  .sonic 
phase  boundaries  and  that  the  accompanying  non-uniqueness  is  an  instability. 

After  dealing  with  various  preliminary  issues  in  Sections  2-4  we  turn,  in  Section  5, 
to  deduce  the  general  solution  forms  to  the  Riemann  problem  that  are  consistent  with 
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the  entropy  inequality.  In  Section  6  we  explicitly  construct  all  solutions  to  the  Hiemann 
problem  corresponding  to  initial  data  with  strains  in  the  low-stain  phase;  the  results  are 
summarized,  and  the  nature  of  the  associated  non-uniqueness  characterized,  in  Section 
7.  Finally  in  Section  8  we  introduce  the  nucleation  criterion  and  the  kinetic  relation, 
thus  leading  to  uniqueness. 

2.  Background.  Consider  longitudinal  motions  of  an  elastic  bar  that  is  regarded 
as  a  one-dimensional  continuum  with  unit  cross-sectional  area.  The  motion  of  the  bar 
is  assumed  to  take  place  isothermally.  During  such  a  motion,  the  particle  at  x  in  the 
reference  configuration  is  carried  to  x  +  u{x,  t)  at  time  t,  where  u{x,  t)  is  the  displacement. 
The  displacement  is  assumed  to  be  continuous  with  piecewise  continuous  first  and  second 
derivatives  throughout  the  regions  of  space-time  to  be  considered.  The  strain  and  particle 
velocity  are  defined  by  7  =  and  v  =  ut  at  points  (x,t)  where  the  derivatives  exist. 
Necessarily,  'y{x,t)  >  —  1  in  order  to  ensure  that  the  mapping  x—*x  u{x,t)  is  invertible 
at  each  instant  t.  The  stress  <7(1,  t)  is  related  to  the  strain  through 

cr  =  <T(7),  (1) 

where  a  is  the  stress  response  function  of  the  material.  At  points  where  7  and  v  are 
smooth,  balance  of  momentum  and  kinematic  compatibility  require  that 

=  pVt,  =  (2) 

where  the  constant  p  is  the  mciss  density  in  the  reference  configuration.  If  there  is  a 
moving  discontinuity  at  x  =  s(t),  the  following  jump  conditions  must  hold: 

a(7)  -  <t(7)  = -pi  (v  -  u),  (n  -  r)  = -i  (7  -  7),  (3) 

where  for  any  function  g{x,t)  we  write  5  =  p(s(t)±,t)  for  the  limiting  values  of  g  on 
either  side  of  the  discontinuity. 

Consider  the  motion  of  the  piece  xi  <  x  <  X2  of  the  bar  during  a  time  interval  [ti,  t2]- 
Suppose  that  7  and  v  are  smooth  on  [xi,X2]  x  [ti,<2]  except  at  the  moving  discontinuity 
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X  =  s(f).  Let  E{t)  be  the  total  mechanical  energy  at  time  t  associated  witli  this  piece  of 
bar; 

Jxi  I 

where  ^(7)  is  the  strain  energy  per  unit  reference  volume  of  the  bar,  i.e. 

^^'"(7)=/  7  >  -1.  (5) 

Jo 

The  following  work-energy  relation  can  be  readily  established; 

cr(x2J)v(x2,t)  -  a(xj,t)v(xiJ)  -  £(t)  =  (6) 

where  the  driving  force  (or  driving  traction)  fit)  acting  on  the  strain  discontinuity  is 
defined  by 

+ 

/  =  fir  7)  =  »il)d-,  -  i(^(7)  +  #(7))(7  -  7),  (7) 

Note  that  /  may  be  interpreted  geometrically  as  the  difference  between  the  area  under 
the  stress-strain  curve  between  7  =7  and  7  =7  and  the  area  of  an  associated  trapezoid 
having  the  same  base.  The  right-hand  side  of  (6)  represents  the  instantaneous  dissipation 
rate  due  to  the  moving  discontinuity  and  the  requirement  that  it  be  non-negative  implies 
that 


/{t)i(0>0.  (8) 

Under  isothermal  conditions,  the  inequality  (8)  is  a  consequence  of  the  second  law  of 
thermodynamics. 

A  motion  of  the  bar  is  governed  by  the  field  equations  (2)  at  all  points  of  smoothness, 
and  the  jump  conditions  (3)  and  the  entropy  inequality  (8)  at  all  discontinuities. 

3.  Material.  Shock  waves  and  phase  boundaries.  In  this  paper,  we  consider  a 
material  whose  stress  response  function  (7(7)  is  twice  continuously  differentiable  with  a 
first  increasing  with  increasing  7,  then  decreasing,  and  finally  increasing  again  as  shown 
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in  Figure  1.  More  specifically  we  suppose  that  there  are  three  nunibers  *viv/,  and  '7,„ 


with 

0  <  7M 

<  7.n  < 

7m  such 

that 

^  >0, 

—  1  <  7 

<  7m, 

=  0, 

7  =  7Afi 

a'(7)  ^ 

<0, 

7m  <  7 

7m , 

i  =0, 

7  =  7  m, 

1  >  0, 

7  >  7m, 

and 

<0, 

—  1  <  7 

7in, 

a"(7)  ■ 

=  0, 

II 

3 

.  >0, 

7  >  7.n 

Moreover,  we  suppose  that  ^(0)  =  0  and  that 
5-(-y)  _oo,  cr'{‘i)  — >  00  as7  — >  —1; 
^(7)  = /^oo7  +  <7r  +  o(l)  as7-->oo, 


(9) 


(10) 


(11) 

(12) 


where  fioo{>  0)  and  ax  are  constants.  The  stress-strain  curve  therefore  consists  of  three 
branches,  two  of  which  are  rising,  while  the  other  is  declining;  it  htis  a  single  inflection 
point  at  the  strain-level  7  =  7i„  and  is  asymptotic,  at  large  tensile  strains,  to  the  straight 
line  a  =  /ioo7  +  <7j.  It  is  useful  for  later  purposes  to  note  that  there  are  two  unique  values 
of  strain  R^o  and  P^o  such  that 

+  or  =  <7(/?oo),  ^'(F’oo)  =  /ioo,  where  -\  <  <  P^  <  7^;  (13) 

Roo  is  the  strain-level  at  which  the  asymptote  a  =  -f  crj  intersects  the  first  branch 
of  the  stress-strain  curve,  while  Poo  is  the  value  of  strain  on  the  first  branch  at  which  the 
slope  equals  the  slope  of  this  asymptote.  Certain  other  material  parameters  are  defined 
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in  the  figure.  In  particular,  the  Maxwell  stress  Oo  is  the  stress-level  for  which  the  two 
hatched  areas  of  Figure  1  are  equal. 

We  shall  say  that  a  particle  of  the  bar  labeled  by  x  in  the  reference  state  is  in  the 
low-strain  phase,  the  ‘"unstable  phase”  or  the  high-strain  phase  at  time  t  during  a  motion 
if  'y{x,t)  lies  in  the  respective  intervals  (  — l,'y,w],(yAf,'ym)  or  [7^,00).  At  a  moving 
discontinuity  i  =  s(^),  the  jump  conditions  (3)  imply 

•  2  ^(7) -a 

Ps  =  -  + 

7-7 


(7) 


(>0). 


(14) 


A  discontinuity  is  called  either  a  shock  wave  or  a  phase  boundary  according  to  whether  7 
and  7  both  lie  in  the  same  phase  or  in  distinct  phases.  The  sound  speed  of  the  material 
at  a  strain  7  is  defined  by 

0(7)  =  /^.  (15) 

where  it  is  necessary  that  7  in  (15)  not  belong  to  the  unstable  phase.  Let  —  -^poolp- 
Let  c  =  c(7)  stand  for  the  sound  speeds  on  the  two  sides  of  a  discontinuity.  The  prop¬ 
agation  speed  i  of  the  discontinuity  is  said  to  be  subsonic  if  [i]  <c  and  c,  intersonic  if 
c<  [i]  <c  or  c<  ji]  <c,  and  supersonic  if  |i|  >c  and  c. 

We  shall  speak  of  a  low-strain  shock  wave  and  a  high-strain  shock  wave  according  to 
whether  the  strains  7, 7  both  belong  to  the  low-strain  phase  or  to  the  high-strain  phase. 
For  the  material  (9)-(12)  considered  here,  it  can  be  readily  seen  that  all  shock  waves  are 
intersonic.  Moreover,  it  follows  from  (7)  and  (9)-(12)  that  the  entropy  inequality  (8) 
holds  at  a  shock  wave  if  and  only  if 


f  7>7 

Low-strain  shock:  <  _ 

if  i  >  0, 

(16) 

[  7<7 

if  i  <  0, 

I  7<7 

High-strain  shock:  <  _ 

if  i  >  0, 

(17) 

[  7>7 

if  s  <  0. 

This  implies  in  particular  that  a  shock  wave  always  moves  into  the  phase  whose  sound 
speed  is  smaller  than  the  speed  |i|  of  the  shock  wave. 
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Turning  next  to  phase  boundaries,  we  will  show  in  the  next  section  that  a  phase 
boundary  for  which  either  7  or  7  is  in  the  unstable  phase  cannot  arise  in  the  problem 
to  be  considered  here.  Thus  suppose  that  7  belongs  to  the  high-strain  phase  and  7  to 
the  low-strain  phase.  In  the  7,  7-plane,  the  set  of  all  pairs  (7,  7)  for  which  7  is  in  the 
high-strain  phase,  7  is  in  the  low-strain  phase  and  the  right  side  of  (14)  is  non-negative 
is  represented  by  the  union  T  of  the  hatched  regions  in  Figure  2;  it  is  the  region  bounded 
by  the  lines  7=  —1,7=  7a/,  7=  7m  and  the  curve  <7(7)  =  <7(7).  The  region  F  of  the 
7,7-plane  will  play  a  major  role  in  the  analysis  that  follows  in  the  next  sections.  The 
boundary  segment  £  is  defined  by 

£:  <7(7)  =  (7(7)  O  7=T(7),  7m<7<7;3, 

where  the  material  parameters  7^3  and  7m  are  defined  in  Figure  1.  By  (14),  i  =  0  at 
points  on  £  and  so  this  segment  represents  instantaneously  stationary  or  equilibrium 
states  of  the  phase  boundary.  One  can  verify  that  T’'(7)  >  0  for  7m  <  7  <  7/3  and  that 
T'(7m)  =  0,  T\^0)  =  cx>;  the  curve  £  therefore  rises  monotonically  as  7  increases.  Next, 
consider  the  curve  T  which  is  defined  as  the  set  of  points  (7,7)  at  which  the  driving  force 
/(7,7)  introduced  in  (7)  vanishes: 

T-.  /(i7)  =  0  ^  7=Q(7),  7  >  703,  (19) 

where  the  material  parameter  703  is  defined  in  Figure  1.  One  can  verify  that  Q'{j)  <  0 
and  that  Qi'y)  —*  Roo,  Q'il)  — >  0  as  7  — >  00  where  Poo  is  the  value  of  strain  defined 
previously  in  (13).  The  curve  T  therefore  declines  monotonically  as  7  increases  as  shown 
in  the  figure.  In  view  of  (8),  a  pheise  boundary  associated  with  a  point  on  T  propagates 
without  dissipation.  Since  /(7,7)  >  0  above  the  entropy  inequality  indicates  that 
i  >  0  there;  likewise  /  <  0  and  i  <  0  below  T.  Consider  next  the  curves  S  and  S  which 
are  defined  as  the  (“sonic”)  curves  on  which  the  speed  h  of  the  phase  boundary  is  equal. 
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respectively,  to  the  sound  speeds  c  and  c: 


+ 

/  +  - 

/  + 

s- 

^(d(7)-d(7))/(7 

~  7)  =  1 

/^'(7) 

7=  P(7),  7  >  70, 

S: 

7((t(7)  -d(7))/(7 

-  7)  =  1 

A'(7) 

7=  P(7),  7  >  7m 

(20) 


One  can  verify  that  ^'(7),  R'i^)  <  0  and  that  P(7)  — >  Poo,  R{l)  Poo,  P^(7  ),  R'il  )  0 

as  7  — »  00  where  Poo  and  Roc  were  defined  earlier  in  (13).  The  curves  $  and  S  therefore 
decline  monotonically  with  increasing  7  as  shown  in  the  figure.  One  can  also  verify  that 
the  three  curves  S,  S  and  do  not  intersect  each  other;  necessarily,  the  curves  P  and 
S  approach  each  other  asymptotically  as  7  — >  00.  The  region  T  is  thus  divided  into  four 
subregions  Ti,  Fa,  Fa  and  F4  by  these  curves.  The  regions  Fi  and  F2  correspond  to  phase 
boundaries  which  propagate  into  the  high-strain  phase  at,  respectively,  intersomc  and 
subsonic  speeds;  the  regions  F3  and  F4  correspond  to  phase  boundaries  which  propagate 
into  the  low-strain  phase  at,  respectively,  subsonic  and  intersonic  speeds.  Points  on  the 
curves  S  and  S  correspond  to  sonic  pha^e  boundaries.  For  the  material  (9)-(12)  under 
consideration,  supersonic  phase  boundaries  cannot  occur.  We  note  that  the  figure  has 
been  drawn  for  the  case  7^  >  Poo  though  we  do  not  assume  this  in  the  analysis. 

Finally  we  consider  the  mapping  (7,7)  — ►  («,/)  defined  by  (7),  (8)  and  (14).  One 

can  verify  that  the  Jacobian  determinant  of  this  mapping  vanishes  when  7,7  corresponds 

to  a  sonic  phase  boundary,  i.e.  on  the  curves  S  and  S-  Considering  the  subsonic  and 

intersonic  regions  separately,  one  can  map  each  of  the  regions  F,  into  the  s, /-plane; 

Figure  3  shows  the  images  F(  that  result  from  this  mapping.  Each  of  the  curves  S\  S\ 

A4'  and  Af'  rises  monotonically  as  s  increases.  (  The  curves  M'  and  A4(.  are  the  images 

+ 

of  the  straight  lines  M  and  M.  shown  in  Figure  2.)  As  i  — ^  Cqo  the  curves  S'  and  M' 
rise  without  bound;  when  i  —>■  ~Coo  the  curve  S'  declines  without  bound.  In  Figure  3, 
Coc  =  yjt^oolp,  Im  =  and  fm  =  /(7m, 7a)-  Note  that  though  this  mapping  of 

F  from  the  7, 7-pIane  to  the  i, /-plane  is  not  one-to-one,  when  restricted  to  the  subsonic 
region  F2  U  F3,  it  is  one-to-one. 

4.  The  Riemann  Problem.  Wave  Fans.  We  now  formulate  the  Riemann  problem 
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for  the  field  equations  and  jump  conditions  (2)-(3).  We  seek  weak  solutions  of  tlie  differ¬ 
ential  equations  (2)  on  the  upper  half  of  the  /-plane  that  satisfy  the  following  initial 
conditions: 


7(a:,0),u(a:,0)  =  < 


—  OQ  <  X  <  0, 


(21) 


[  0  <  z  <  oo, 

where  and  vr  are  given  constants  with  7l  >  —1  and  7/?  >  1. 

Since  the  initial  value  problem  described  above  is  invariant  under  the  scale  change 
t—^kt,  x^kx,  we  restrict  attention  to  solutions  that  have  this  property  as  well.  Where 
such  solutions  exist,  they  must  have  the  form  7(x, /)  ==  —  ^’(0  where  ^  =  x/t. 

It  then  follows  from  (2)  that  for  such  solutions,  either  7  and  v  are  both  constant,  or  they 
are  the  wave  fans  which  are  given  by 


(22) 


v'iO  =  -^7(0- 


(23) 


The  left  side  of  (22)  must  necessarily  be  non-negative  for  any  7  that  satisfies  (22).  Thus 
for  the  material  (9)-(12)  wave  fans  can  occur  only  if  7  takes  values  in  either  the  low- 
strain  phase  or  the  high-strain  phase.  We  shall  speak  of  a  low-strain  fan  or  a  high- 
strain  fan  according  to  whether  7  belongs  to  the  low-strain  phase  or  the  high-strain 
phase  respectively.  In  view  of  (9)-(12),  it  follows  that  (22)  defines  a  unique  function 
7(0  =  7^^^(0  for  -00  <  (  <  00,  such  that  7<^^  €  (-1,7m];  7*^^  describes  the  strain  field 
in  a  low-strain  fan.  Similarly  (22)  can  be  uniquely  solved  for  a  function  7(^)  =  7*^^('f) 
for  -Coo  <  i  <  Coo,  such  that  7'^)  €  [7m, 00);  7^^^  describes  the  strain  field  in  a  high- 
strain  fan.  Thus  (22),  (23)  lead  to  the  two  wave  fans 

7(0  =  7'’ne),  i  =  L,H,  where  (24) 

-{v^'\()-v^'\^o)]  =  ±  c{^)d'y,  1  =  L,H,  (25) 

and  ^0  describes  an  arbitrary  ray  xjt  -  io  within  the  fan;  necessarily  must  lie  in  the 
interval  (—00,00)  for  a  low-strain  fan  and  in  the  interval  (— Coo,Coo)  for  a  high-strain 
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fan.  Ttie  positive  and  negative  signs  are  taken  in  the  right  side  of  (25)  according  to 
whether  the  wave  fan  occurs  in  the  first  quadrant  or  second  quadrant  of  tiie  j-.t-piain*, 
respectively. 

Consider  two  rays  j  =  ct  and  x  =  <•/  with  c  <  c  in  the  (j'.f)-plane  hetwecui  wliich 
the  field  is  a  fan.  Let  5=  +  t’=  v{ct  +  J),  ~/=  r~  v{c(-.t)  denote  the 

limiting  values  from  within  the  fan  of  strain  and  particle  velocity  at  these  rays.  It  follows 
from  (15)  and  (22)  that 

c  —  ±c{l),  c  =  ±c(5),  (20) 


and  from  (23)  that 


r 


±  I  c(7)<f7. 
J  > 


(27) 


where  the  positive  and  negative  signs  are  taken  according  to  whether  the  fan  occurs  in 
the  first  or  second  quadrant  of  the  x,f-plane.  respectively.  Equation  (27)  is  the  analog 
for  a  fan  of  the  kinematic  jump  condition  for  a  discontinuity  in  (3).  Since  the  field  within 
the  fan  is  smooth,  the  entropy  inequality  is  trivially  satisfied  at  points  within  it.  For  the 
material  (9)-(12),  it  follows  from  (22)  that  necessarily 

Low-strain  fan: 

High-strain  fan; 


4-  “ 

7<7 
+  - 

if  fan  is  in  first  quadrant. 

(28) 

7  >5 

if  fan  is  in  second  quadrant. 

4*  — 

7>7 

if  fan  is  in  first  quadrant. 

(29) 

—<>4' 
A 
->  1 

if  fan  is  in  second  quadrant. 

Equations  (28)-(29)  are  the  analog  for  fans  of  equations  (16)-(17)  for  shotAS.  Conversely, 
given  numbers  (7,t;),  (7,n)  which  conform  to  (27)-(29).  one  can  construct  a  unique  fan 
between  the  rays  x  ~  ct  and  x  —  ct  where  c  is  given  by  (26). 

The  general  scale-invariant  solution  to  the  Riemann  problem  has  the  form  shown  in 
Figure  4;  between  any  two  rays  x  =  s,t  and  x  =  s.+ii  the  fields  7,  c’  are  either  constants 
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or  fans;  the  rays  themselves  may  or  may  not  correspond  to  discontinuities.  If  x  =  s,/  is 
a  discontinuity,  the  jump  conditions  (3)  must  be  satisfied  across  it  so  that 

{s,(7,  -  7,_i)  =  -(u.  - 

(30) 

pi?  =  (<t(7.)  -  5-(7._i))/(7.-  -  7,_i), 

where  {'yi,Vi)  and  (7,_i,Ui_i)  are  the  limiting  values  of  strain  and  particle  velocity  on  the 
right  and  left  respectively  of  this  discontinuity.  Let  /,  =  /(7,-i,7,)  with  /  defined  by  (7) 
stand  for  the  driving  force  on  this  discontinuity;  the  entropy  inequality  (8)  then  requires 
that 


/.i,  >0.  (31) 

An  admissible  solution  of  the  Riemann  problem  is  a  pair  '){x,t),v{x,t)  of  the  form  just 
described  with  (30)-(31)  enforced  at  all  discontinuities. 

5.  The  Structure  of  Admissible  Solutions  to  the  Riemann  Problem.  We  shall 
say  that  the  initial  data  in  (21)  is  metastable  if  neither  of  the  initial  strains  belong 

to  the  unstable  phase.  Before  constructing  explicit  global  solutions  to  the  Riemann 
problem,  it  is  helpful  to  establish  some  general  results  pertaining  to  the  permissible 
solucion  forms  that  are  consistent  with  the  entropy  inequality. 

Let  (7,u)  be  an  admissible  solution  of  the  Riemann  problem  with  metastable  initial 
data. 


(i)  The  strain  7(1,  t)  does  not  belong  to  the  unstable  phase  at  any  point  (x,t) 
in  the  upper-half  plane. 

This  result  implies  that  if  the  initial  data  does  not  involve  the  unstable  phase,  then 
at  no  later  time  does  the  solution  involve  the  unstable  phase.  We  prove  this  claim 
by  contradiction.  Suppose  that  this  proposition  is  false.  Since  unstable  phase  fans 
and  shocks  do  not  exist,  there  must  necessarily  be  two  phase  boundaries  in  the  {x,t)- 
plane,  say  x  =  Skt  and  x  —  with  i*  <  ijt+i,  such  that  the  state  between  them 

is  constant  with  the  associated  strain  7*  in  the  unstable  phase  and  with  neither  of  the 
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strains  7^.1  =  'y{skt  —  ,t)  and  nik+i  =  l{sk+ii+,t)  in  the  unstable  phase.  From  (7)  and 
(9)-(12)  one  sees  that  the  driving  force  on  the  discontinuity  x  =  Skt  is  positive;  the 
entropy  inequality  (31)  thus  implies  that  Sk  >  0.  Similarly  the  driving  force  on  i  =  s^+it 
must  be  negative  and  so  i^+i  <  0-  Thus  Sk  >  i/t+i  which  is  a  contradiction.  This 
establishes  the  proposition. 

The  next  six  propositions  are  concerned  with  the  possibility  of  having  shock  waves, 
phase  boundaries  and  wave  fans  adjacent  to  each  other.  When  addressing  fans,  it  is 
suflBcient  for  our  purposes  to  restrict  attention  to  fans  which  do  not  terminate  on  discon¬ 
tinuities.  Thus  in  the  propositions  that  follow,  if  x  =  Skt  and  x  =  Sk+it  are  two  rays  in 
the  (x,f)-plane  between  which  the  field  is  a  fan,  these  two  rays  will  be  assumed  not  to 
be  discontinuities.  Features  of  fans  which  do  terminate  at  discontinuities  can  be  deduced 
from  the  following  results  by  suitable  limiting  arguments. 

Let  (7,u)  be  an  admissible  solution  of  the  Riemann  problem  with  metastable  initial 
data* 


(ii)  Let  X  =  Skt,  x  —  Sk+\t,  x  =  ifc+2t  and  x  =  Sk^zt,  be  four  rays  in  the  same 
quadrant  of  the  (x,  t)-plan€.  If  the  field  in  the  interior  wedge  <  x  <  $k+2i 
is  constant,  then  the  field  in  the  other  two  wedges  cannot  both  be  fans. 

(Hi)  Let  X  =  Skt  and  x  =  be  two  rays  in  the  same  quadrant  of  the  {x,t)- 
plane  between  which  the  field  is  constant.  Then  these  two  rays  cannot  both  be 
shock  waves. 

(iv)  Let  X  =  Skt  and  x  =  Sk+it  be  two  rays  in  the  same  quadrant  of  the  (x,t)- 
plane  between  which  the  field  is  constant.  Then  these  two  rays  cannot  both  be 
phase  boundaries. 

(v)  Let  X  =  Skt,  X  =  and  x  —  Sk+zl  be  three  rays  in  the  same  quadrant 

of  the  {x,t)-plane.  Suppose  that  the  field  between  any  two  of  these  rays  is  a 
fan.  Then  the  third  ray  cannot  be  a  shock. 

(vi)  Let  X  =  Skt,  X  =  i/fc+it  and  x  =  Sk+2t  be  three  rays  lu  the  same  quadrant 
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of  the  {x,t)-plane.  Suppose  that  the  field  between  the  two  slowest  rays  is  a 
fan.  Then  the  remaining  ray  cannot  be  a  phase  boundary.  (The  converse  case 
is  possible:  if  the  field  between  the  two  fastest  rays  is  a  fan  the  remaining  ray 
may  be  a  phase  boundary.) 

(vii)  Let  X  =  Skt  and  x  =  Sk+it  be  two  rays  in  the  same  quadrant  of  the 
{x,t)-plane.  If  the  slower  ray  is  a  shock,  then  the  faster  ray  cannot  be  a  phase 
boundary.  (The  converse  case  is  possible:  if  the  faster  ray  is  a  shock,  the 
slower  ray  may  be  a  phase  boundary.) 

The  first  three  of  these  propositions  state  that  two  wave  fans,  two  shock  waves  and 
two  phase  boundaries  cannot  be  adjacent  to  each  other.  The  next  one  states  that  a 
fan  and  a  shock  cannot  be  adjacent  to  each  other.  On  the  other  hand  according  to 
proposition  (vi),  a  fan  and  a  phase  boundary  may  be  adjacent  to  each  other  provided  the 
phase  boundary  is  subsonic.  Similarly  a  shock  and  a  phase  boundary  may  be  adjacent 
to  each  other  provided  the  phaise  boundary  travels  more  slowly  than  the  shock. 

It  is  clearly  sufficient  to  prove  these  results  in  any  one  quadrant  of  the  upper-half  of 
the  (i,  t)-plane  and  so  we  shall  consider  only  the  first  quadrant.  Thus  in  each  of  these 
propositions  we  have  0  <  i*,  <  ijt+i  <  -Sfc+s  <  ^*+3. 

To  prove  proposition  (ii),  suppose  that  it  is  false  so  that  the  field  in  Skt  <  x  <  Sk+it 
and  ifc+2^  <  x  <  Sk+st  are  fans  while  the  field  is  constant  in  the  intermediate  wedge.  It 
is  not  possible  that  one  of  these  fans  is  a  low-strain  fan  while  the  other  is  a  high-strain 
fan,  since  such  fans  belong  to  different  phases  and  so  must  necessarily  be  separated  by 
a  phase  boundary.  For  two  fans  of  the  same  type,  (22)  directly  shows  that  they  must  in 
fact  be  smoothly  connected  to  each  other  to  become  a  single  fan  with  Sk+i  =  Sk+2-  This 
contradict.,  the  assumption  that  Sk+i  <  ifc+2-  The  assertion  (ii)  is  thus  proved. 

We  now  turn  to  the  proof  of  proposition  (iii).  Suppose  that  the  proposition  is  false 
so  that  X  —  Skt  and  x  =  Sk+it  are  both  shocks  and  the  field  between  them  is  constant. 
Note  first  that  a  low-strain  shock  cannot  be  adjacent  to  a  high-strain  shock,  since  they 
must  be  separated  by  a  phase  boundary.  Suppose  that  both  x  =  Skt  and  x  —  Sk+it  are 
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low-strain  shocks.  The  strain  and  velocity  field  thus  have  the  form 


Ik-liVk-l, 


X  =  i*/-. 


7{i,t),n(x,0  = 


lk,Vk, 


Skt  <  I  < 


(32) 


I  7fc+i)^fc+i?  X  —  Sk+\t-\-, 

where  7^-1,  7fc  and  7^+1  are  distinct  and  all  three  belong  to  (  — l,7Af].  According  to  (31 ), 
the  driving  forces  acting  on  these  two  shock  waves  must  be  non-negative.  'Vith  the  help 
of  (7)  one  finds  that  this  implies  that  7^-1  <  7fc  <  lk+\-  Since  d-'(7)  >  0  and  d"(7)  <  0 
on  (  — 1,7m)i  this,  together  with  the  fact  that  —1  <  ')k-\  <  7*  <  Ik+i  <  Ixt-  implies  that 

^hk+i)  -  ^{ik)  ^  <7(70  -  g(7fc-i)  (33^ 

I'fe-t-i  -  Ik  Ik-  7fc-i 

Equations  (30)  and  (33)  then  yield  i*  >  i/t+i  which  is  a  contradiction.  In  a  similar  way. 
we  can  prove  that  two  high-strain  shocks  cannot  be  adjacent  to  each  other.  This  proves 
proposition  (iii). 

Propositions  (iv)  and  (vii)  may  be  established  by  arguments  that  are  very  similar  to 
the  above. 


Next  we  turn  to  the  proof  of  proposition  (v).  Suppose  that  it  is  fcdse.  Note  first  that 
a  low-strain  shock  cannot  be  adjacent  to  a  high-strain  fan  (or  vice  versa)  since  they  must 
be  separated  by  a  phase  boundary.  Suppose  that  x  —  Skt  is  a  low-strain  shock  and  that 
the  field  between  the  rays  x  =  Sk+ii  and  x  =  sjt+jt  is  a  low-strain  fan: 

{  7fc-i,Ufc-i,  x^Skt-, 


'y{x,t),v{x,t)  =  < 


lk,Vk, 


Skt  <!  X  ^  'Sjk+lt, 


Sk+lt  ^  X  ^  ^k-k^tj 


(34) 


[  7fc+2,Ufc+2.  X  =  +  , 

where  7*-i,7fc,  and  7fc4.2  all  lie  in  the  low-strain  phase  and  7^^^  is  the  strain  field  in 
a  low-strain  fan  given  by  (24).  For  the  material  (9)-(12),  equation  (22)  and  the  fact 
that  ifc+i  <  ifc+2  implies  that  7*  >  7/t+2-  Next,  the  entropy  inequality  (31)  requires 
the  driving  force  f{‘^k~^■,lk)  at  x  =  Skt  to  be  non-negative;  by  (7),  this  implies  that 
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7fc-i  <  Ik-  Since  ^'(l)  >  0  <^"(7)  <  0  011  [7it-iOfc]i  t,his,  together  with  the  fact  that 

—  1  <  7jb-i  <  7a:  <  7m  implies  that 


^'(7-t)  < 


^(7fc)  -<7(7fc-i) 
Ik  ~  Ik-i 


Equations  (15),  (30)  now  yield  Sk  >  i^+i  which  is  a  contradiction. 

The  remaining  cases  where  x  —  ijt+2^  is  a  low-strain  shock  and  the  field  between  the 
rays  x  =  Skt  and  x  =  is  a  low-strain  fan,  and  when  the  shock  ,'nd  fan  are  both 

high-strain  ones,  can  be  treated  similarly.  This  establishes  proposition  (v). 

The  proof  on  proposition  (vi)  is  entirely  analogous. 

The  preceding  results  imply  that  the  form  of  admissible  solutions  to  the  Riemann 
problem  with  metastable  initial  data  is  in  fact  much  simpler  than  that  described  in 
Figure  4.  Note  that  the  results  (iii)-(vii)  depend  criticedly  on  the  entropy  inequality  (31 ). 

6.  Explicit  Solutions  to  the  Riemann  Problem.  The  results  established  in  the 
preceding  section  allow  one  to  determine  all  admissible  solutions  to  the  Riemann  problem 
in  the  case  of  metastable  initial  data.  From  here  on  we  shall  consider  only  the  special 
Riemann  problem  in  which  the  initial  strains  ■yi  and  yn  are  both  in  the  low-strain  phase, 
and  7/j  is  smaller  than 


7l€(-1,7m],  7fl^(-T7M],  7fl  <  7t-  (36) 

At  the  initial  instant,  the  entire  bar  is  in  the  low-strain  phase.  At  a  later  instant,  a 
particle  of  the  bar  may  or  may  not  change  its  phase.  It  is  convenient  in  the  following 
analysis  to  consider  these  two  cases  separately. 

(i)  Solutions  involving  no  phase  change.  In  this  case,  the  solution  does  not  involve 
any  phase  boundaries.  In  view  of  the  first  proposition  in  Section  5,  neither  does  it  involve 
the  unstable  phase  at  any  time  it  >  0.  Next,  in  view  of  propositions  (ii),  (iii)  and  (v),  the 
solution  7,  V  can  only  involve  a  single  low-strain  shock  wave  or  a  single  low-strain  fan  in 
each  quadrant  of  the  upper  half  of  the  x,<-plane.  Thus  the  solution  must  have  one  of  the 
four  forms  shown  in  Figure  5(a)-(d). 
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(i(o.))  Solution  with  two  shock  waves.  Figure  5(a).  Consider  a  solution  having  the 
form  shown  in  Figure  5(a); 


IlA'L, 

—00  <  X  <  Sit, 

Sit  <  X  <  S2t, 

(37 

.  7R,r'ft, 

s^i  <  X  <  00, 

in  which  and  are  to  be  found  such  that  7  €  (-l,7,v/]  and  ii  <  0  <  S2. 

The  jump  conditions  (30)  and  the  entropy  inequality  (31)  at  each  of  the  two  shock 
waves  require  that 


-(uft-J))  =  62(7/? -7), 

.  ^(7ft)-^(7)  ^  - 

<*2  =  _  ,  7h  >  7, 

(38) 

N  pilR-'y) 

-(u-ul)  =  ii(7-7L), 

.  _  <T(7)-^(7z.)  -  , 

•Sl  -  _  ^  ^  <  IL- 

(39) 

\  Pi't-ll) 

The  inequalities  in  (36),  (38),  (39),  together  with  the  requirement  that  7  be  in  the  low- 
strain  phase,  imply  that 

-1  <  7  <  7ft.  (40) 

Combining  (38)-(39)  yields 

= /f(7),  (41) 

where  H{'y)  is  defined  on  (  —  1,7ft]  by 

^(7)  =  -■ v/[<7(7h)  -  ^(7)](7«  -  l)Jp  -  0^(7)  -  ^(7l-)](7  -  7l)/p-  (42) 

It  can  be  verified  that  Hiy)  increases  monotonically  on  (  — 1,7r]  from  the  value  — oc 
at  7  =  -1  to  the  value  Hr  =  if (7/?)  at  7  =  7/?;  see  Lin  [19].  Thus  if  the  initial 
data  is  such  that  —00  <  vr  —  vl  <  Hr,  there  is  a  unique  root  7  of  (41)  in  the  range 
-I  <  7  <  7ft(<  1m)-  The  remaining  unknowns  v,ii  and  ^2  are  then  given  immediately 
by  (38)  and  (39). 
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Thus,  thtrt  exists  a  unique  admissible  solution  of  the  form  (37)  corresponding  to 
Figure  5(a)  if  and  only  if  the  given  initial  data  (21),  (36)  is  such  that  — oc  <  vn  —  vi  < 
Hr. 


(i(b))  Solution  with  a  shock  wave  and  a  wave  fan.  Figure  5(b).  Consider  next  a 
solution  of  the  form  shown  in  Figure  5(b): 

f  -oc  <  X  <  st. 


- 


st  <  X  <  ct. 


ct  <  X  <  crI, 


(  1r,vr, 


CRt  <  X  <  oc, 


(43) 


where  7,  v  s,  c  and  cr  are  to  be  determined  such  that  7  G  (  —  1,7a/)  and  s  <  0  <  c  <  cr. 
The  functions  7(x/t)  and  v{x/t)  are  the  strain  and  velocity  fields  pertaining  to  a  low- 
strain  fan  and  are  given  by  (24)-(25). 

At  the  shock  wave  x  =  ii,  the  jump  conditions  (30)  and  the  entropy  inequality  (31) 
must  hold: 


n  =  ui  -  i(7  -  7Z,),  s  - 


(t(7)  -  a{n) 


7  <  7/-- 


\  />(7  -  7l) 

Turning  next  to  the  fan,  and  by  using  (26)-(28),  one  finds 


(44) 


-{v-vr)=I  c{-f)d-y,  cr  =  c(7r),  c  =  c(7),  7> 

•'T'R 

Equations  (44)-(45)  may  now  be  combined  to  yield 


IR- 


(45) 


vr-vl  =  H[1) 


(46) 


where 

/f(7)  =  -  -  ^(7)](7L  -  7)/^  +  /  c(e)de,  7/?  <  7  <  IL-  (47) 

One  can  verify  that  Ff(7)  increeuses  monotonically  on  [7r,  ')if[  from  the  value  Hr  =  H{^r) 
at  7  =  7R  to  the  value  Hi  =  /f(7L)  at  7  =  7/,;  see  [19].  Thus  if  the  iidtial  data  (21), 
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(36)  is  such  that  Hfi  <  Vft—  vl  <  Hi,  equation  (46)  yields  a  unique  root  1  in  the  interval 
')R  <1  <  71.  The  remaining  unknowns  i\  s,  c  and  cr  are  then  given  by  (44),  (45). 

Thus,  there  is  a  unique  admissible  solution  of  the  form  (43)  corresponding  to  Figure 
5(b)  if  and  only  if  the  given  initial  data  (21),  (36)  is  such  that  Hr  <  i  -R  -  VL  <  Hl. 

(i(c))  Solution  with  two  wave  fans.  Figure  5(c).  We  seek  solutions  in  the  form  of 
Figure  5(c): 


—  OC  <  X  <  CLt, 

7](x/f),Vi(x/0, 

Clt  <  X  <  c't. 

c't  <  X  <  Ct, 

Ct  <  X  <  CrI, 

iRi Vr, 

crI  <  X  <  00, 

(48) 


where  7,  v,  ci,  c',  c  and  cr  are  to  be  determined  such  that  7  €  (— l,7Af]  and  c^  < 
c'<  0  <  c  <  cjcj.  The  functions  71,  Ui,  72,  V2,  are  the  fields  that  are  appropriate  to  a 
low-strain  fan  and  are  given  by  (24)-(25). 

The  analysis  of  this  Ccise  is  similar  to  that  of  the  preceding  one.  One  finds  that  there 
is  a  unique  admissible  solution  of  the  form  (48)  corresponding  to  Figure  5(c)  if  and  only 
if  the  given  initial  data  (21),  (36)  is  such  that  Hl  <  vr  —  vl  <  H\i  where 

H{^)=  f  c{^)d-f  +  I  c{-t)d'f,  (49) 

J-IR  JlL 

and//i  =  //(7L),/fAf  =  /^(7m). 

(i(d))  Solution  with  a  shock  wave  and  a  wave  fan.  Figure  5(d).  For  initial  strains  in 
the  low-strain  phase  with  7n  <  7l  as  considered  here,  one  finds  that  solutions  having 
the  form  of  Figure  5(d)  do  not  exist.  In  the  reverse  case  7/?  >  7^,,  one  finds  that  such 
solutions  do  exist  in  place  of  solutions  of  the  form  in  Figure  5(b)  which  now  do  not  exist. 

It  is  readily  seen  from  (42), (47)  and  (49)  that  f/(7)  is  continuous  at  7  =  7/?  and 
7  =  In  the  respective  limits  vr  —  vi  Hr—,  and  vr  —  vl  —*  Hr+,  the  solution 
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forms  shown  in  Figure  5(a)  and  Figure  5(b)  coincide.  The  limiting  solution  involves  a 
single  shock  wave  traveling  left  and  has  no  other  shock  waves.  Similarly,  the  solution 
forms  :n  Figure  5(b)  and  Figure  5(c)  coincide  in  the  respective  limits  Vfi  —  =  //l-. 

t’H  —  vi  —  Hi+.  In  this  case  the  solution  involves  a  single  rightward  moving  wave  fan. 

In  summary,  when  initial  data  (21),  (36)  is  such  that  — oc  <  vn  —  vi  <  Hm.  it  has 
been  shown  that  there  is  a  unique  admissible  solution  that  involves  no  phase  change  to 
the  Riemann  problem.  Further  discussion  of  these  solutions  is  postponed  until  subsection 
(iii). 

In  order  to  find  solutions  when  the  initial  data  is  such  that  vn  —  vi  >  Hm  we  must 
consider  solutions  which  involve  a  phase  change. 

(ii)  Solutions  involving  a  phase  change.  If  a  solution  involves  a  phase  change,  the 
results  (iv),  (vi)  and  (vii)  of  Section  5  show  that  each  quadrant  of  the  upper-half  of  the 
(x,  t)-plane  has  precisely  one  phase  boundary  together  with  either  one  shock  wave  or  one 
wave  fan.  Moreover,  each  phase  boundary  is  necessarily  subsonic  so  that  the  speed  of  the 
shock  wave  or  wave  fan  is  greater  than  that  of  the  phase  boundary.  Thus  the  solution 
must  have  one  of  the  four  forms  shown  in  Figure  6(a)-(d). 

(ii(a))  Solutions  with  two  phase  boundaries  and  two  fans.  Figure  6(a).  Suppose  that 
the  solution  has  the  form  shown  in  Figure  6(a): 


1L,VL, 

—  00  <  X  <  Clt, 

7i(x/<),t)i(x/<), 

Clt  <  X  <  c^t. 

7a, va. 

CaI  <  X  <  Sit, 

7,u, 

Sit  <  X  <  S2t, 

+  + 

7,u, 

S2t  <  X  <  %t. 

72{xlt),V2{x(t), 

ct  <  X  <  Cfit, 

7r,  Vr, 

Cjit  <  X  <  00, 

where  7>t,Uy4i  ci,  c^,  c,  C/j, ii  and  S2  are  to  be  determined  such  that  74,7  € 
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(  — I'.TA/]:  ^  [7m, oo)  and  c/,  <  c>i  <  ii  <  0  <  i-i  <  <^  <  (^R-  The  functions 

7i(x/f), {)i{x/<),72(^/0i ^’2(-^'/0  fields  appropriate  to  a  low-strain  fan  and  are 

given  by  (‘24)-(25). 

The  jump  conditions  (30)  at  the  two  phase  boundaries  lead  to 


V  =  V  —  S2{if  ~  7),  ^2  ~ 


^(7)  -a(7) 
^  M7-7) 


v  =  -  ii(7  -  7^),  -Si  = -, 


<t(7)  -^(7.4) 


N  pO-1a) 

while  the  requirements  (26),  (27)  at  the  two  fans  give 


(51) 


(52) 


+ 

r7 


^  -f*  +  f 

c  -  c(7),  Cfl  =  c(7fl),  ~(v  -  Ufl)  =  /  c(7)(i7, 

JyR 

pA 

ca  =  -c(7>i),  Cl  =  -c(7l),  va-vl=  /  €(7)^7. 


(53) 

(54) 


The  roots  of  these  equations  are  subject  to  the  restrictions  imposed  by  the  entropy  in¬ 
equality  (31)  at  each  phase  boundary,  the  inequality  (28)  at  each  fan,  and  the  requirement 
that  the  strains  belong  to  either  low-strain  phase  or  high-strain  phase  as  appropriate. 
These  restrictions  lead  to  the  following  inequalities 


/(7,7)>0,  /(7^,7)<0, 


7  >  n/R,  lA  >  7L, 


-1  <  7  <  7Af, 


T  <  7a  <  7m,  7  >  7,n- 


(55) 

(56) 

(57) 


The  10  equations  (51)-(54)  are  to  be  solved  for  the  12  unknown  quantities  listed  below 
(50).  Therefore,  one  anticipates  that  when  there  exists  a  solution,  there  would  in  fact  be 
a  two-parameter  family  of  solutions]  see  James  [12].  For  reasons  of  algebraic  simplicity, 
we  assume  that 

S2  =  -Si  =  i,  (58) 
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where  i  >  0.  This  assumption  does  not  change  the  essential  characteristics  of  the  so¬ 
lutions,  but  leads  to  a  considerable  simplification  in  the  analysis.  In  fact,  it  is  possible 
to  show  (as  in  Abeyaratne  and  Knowles  [3])  that  (58)  necessraily  holds  once  a  kinetic 
relation  is  imposed.  One  now  expects  a  one-,  rather  than  a  two-,  parameter  family  of  so¬ 
lutions.  In  what  follows,  the  strain  7  snail  be  treated  as  this  parameter.  The  assumption 
(58)  also  leads  to  7^1  =7. 

It  follows  from  (55)-(58)  and  the  requirement  0  <  s  <c  in  (50)  that 

c(7)>6,  /(7,7)>0,  7>7l,  d(7)>(T(7)  (59) 

necessarily  hold.  These  inequalities  define  a  region  D\  C  Fa  in  the  7,7-plane;  Figure  7 
displays  this  region  in  the  case 

7ot  ^  T’oq  >  7L)  7h  ^  Root 

and  for  reasons  of  definiteness  we  shall  frame  our  analysis  from  hereon  for  this  particular 
case;  the  analysis  can  be  trivially  modified  to  handle  the  case  when  the  ordering  of  the 
strains  differs  from  (60). 

By  combining  equations  (51)-(54),  one  can  reduce  the  question  of  their  solvability  to 
the  following  problem:  given  7,  find  a  root  7  with  (7,7)  €  Di  of  the  equation 

Ufi-UL  =  G(7,7),  (61) 

where  G  is  defined  for  (7,7)  €  D\  by 

+  +  - - 

G(7,7)s  Fc(7)<i7+  /  c(7)</7  +  2VW7)-a(7)l(7-7)//>-  (62) 

J'lR 

The  curves  AD,  AI  and  CD  of  Figure  7  which  comprise  part  of  the  boundary  of  the 
region  Di  were  defined  previously  in  (18)-(20).  At  each  fixed  7  >  703,  the  function  G 
in  (62)  increases  monotonically  with  increasing  7.  Therefore,  for  each  7  >  703,  (61)  can 
be  solved  uniquely  for  7  provided  vr  —  lies  in  a  suitable  range;  for  703  S  '^  <  7/?, 
this  range  is  G(7,(5(7))  <  vr  —  vi  <  G{l,T{‘f))]  for  7/?  <  7  <  7/,  it  is  G{'y,Q{l))  < 
VR-  vl  <  G(7,P(7));  and  for  7  >  7/,  it  is  ^7(7,71.)  <  vr  -  vl  <  G{i,P{l)).  Here 
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li{>  Im)  is  the  strain-level  at  which  the  driving  force  /(7/,7z,)  vanishes,  see  Figure  7. 
Once  7  has  been  thus  determined,  all  of  the  other  unknowns  can  be  found  in  terms  of 
7  and  the  initial  data  without  further  restriction.  Further  discussion  of  this  solution  is 
postponed  until  subsection  (iii). 

(ii(b))  Solutions  with  two  pha^e  boundaries,  one  shock  wave  and  one  wave  fan.  Figure 
6(b).  Consider  next  solutions  having  the  form  of  Figure  6(b): 

-OC  <  X  <  Sit, 

IAiVa,  Sit<X<S2t, 


7,w, 

S2t  <  X  <  s^i. 

-(■  -(■ 

7,u, 

•S3t  <  ^  <  Cl, 

l,v, 

ct  <  X  <  CRt, 

1R,vr, 

CRt  <  X  <  00, 

(63) 


where  'yAiVA,  7,  v,  7,  v,  si,  S2,  S3,  c  and  cr  are  to  be  determined  such  that  7^,  7  6  (— 1, 7;V/], 
7  €  [7m,  00)  and  ii  <  i2  <  0  <  is  <  c  <  ch-  The  functions  7,  v  correspond  to  the  fields 
in  a  low-s'rain  fan  and  are  given  by  (24)-(25). 

The  analysis  of  this  case  is  similar  to  the  previous  one  and  so  we  merely  state  the 
results.  Again,  there  is  a  two-parameter  family  of  solutions  in  general,  and  for  algebraic 
simplicity  we  assume  that  33  =  — =  s,  thus  reducing  it  to  a  one-parameter  family  of 
solutions. 

The  strains  (7,7)  must  lie  in  the  region  D2  C  Fs  shown  in  Figure  7.  Given  the  strain 
_  ^ 

7,  one  is  to  find  a  root  7  with  (7,7)  €  D2  of  the  equation 

VR-VL  =  Gii,l),  (64) 


where  G  is  defined  by 


G(7,7)  =  -\/[d-(7i)  -  d(7)](7L  -  7)/p  +  /  €(7)^7 
+2J [0(7)  -  d(7)](7  -  7)/p,  for  (7, 7)  € 


(65) 


23 


—  -I- 

At  each  fixed  7  >  7/  the  function  G  increases  inonotonically  with  7.  1  herefore  for 
each  7  in  the  respective  ranges  7/  <  7  <  7a-  and  7  >  7a,  (64)  can  be  uniquely  solved 
for  7  provided  G{l,Q{l))  <  vr  -  Vi  <  G(7,7i,)  and  G(7,7a)  <  vr  -  vi  <  6'(7,7a) 
respectively.  Here  7a'(>  7m)  is  the  value  of  strain  at  which  the  driving  force  /(7A',7/t) 
vanishes,  see  Figure  7.  Once  7  has  been  thus  determined,  all  of  the  other  unknowns  can 
be  found  in  terms  of  7  and  the  given  initial  data  w-ithout  further  restriction.  Further 
discussion  of  this  solution  is  postponed  until  subsection  (iii). 

(ii(c))  Solutions  with  two  phase  boundaries  and  two  shock  waves.  Figure  fi(c).  Here 
we  seek  solutions  in  the  form  of  Figure  6(c); 


—  00  <  X  <  Sit, 

Sit  <  X  <  s^t. 

7,^, 

S2t  <  X  <  Sjt, 

+  + 

Szt  <  X  <  54^, 

7,^^, 

.  IR.vr, 

s^t  <  X  <  00, 

(66) 


in  which  7a,  7,  u,  7,  n,  51,62,53,34  are  to  be  determined  such  that  7^,  7  €  (  — l,7Af], 

7  G  [7m,  00)  and  ii  <  ^2  <  0  <  53  <  54. 

The  analysis  of  this  case  is  again  similar  to  the  previous  ones.  There  is  again  a  two- 

parameter  family  of  solutions  in  general,  and  for  algebraic  simplicity  we  assume  that 

53  =  —52  =  5,  thuvS  reducing  it  to  a  one-parameter  family  of  solutions. 

The  strains  (7,7)  must  lie  in  the  region  D3  C  Fa  shown  in  Figure  7.  Given  the  strain 
-  +  -  + 

7,  one  is  to  find  a  root  7  with  (7,7)  €  D3  of  the  equation 

vr-vl  =  G{i,^)  (67) 


where  G  is  defined 
G(7,7)  =  -  V 


by 

'^['^(7£-)  -  <5-(7)](7l  -  7)/p  -  \j{cr{lR)  -  ~  ^)/P 

-h2y[<T(7)-d(7)](7-7)/>,  /or(7,7)  €  £>3. 


(68) 
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At  each  fixed  1  the  function  G  increases  nionotonically  with  increasing  I  and 

therefore  (67)  can  be  uniquely  solved  for  1  provided  Gil.Qi  7))  <  VH  -  Cl  <  G{1.',h}. 
Once  7  has  been  determined,  all  of  the  other  unknowns  can  be  found  without  further 
restriction. 

(ii(d))  Solutions  with  two  phase  boundaries,  one  fan  and  one  shock  wains.  Ftguri 
6(d).  One  can  show  that  solutions  having  the  form  shown  in  Figure  6(d)  do  not  exist  for 
the  initial  data  (21),  (36)  being  studied  here.  Such  solutions  do  exist  in  the  case  -jn  > 
in  which  event  solutions  of  the  form  of  Figure  G(b)  do  not  exist. 

It  is  clear  from  (62).  (65)  and  (68)  that  G  is  continuous  across  the  lines  IJ  and  K L 
in  the  (7,7)-plane  shown  in  Figure  7.  In  the  respective  limits,  for  >  7/.  vh  —  i'l  —* 
G{1,')l~)  and  Vr  —  vi  — +  6''(7,7c  +  )<  the  solution  forms  shown  in  F'igure  6(a)  and  Figure 
6(b)  coincide.  This  limiting  solution  involves  a  wave  fan  traveling  right  and  two  pha.se 
boundaries.  Similarly,  the  solution  forms  in  Figure  6(b)  and  Figure  6(c)  coincide  in  the 
limits  —  — »  G(7,7i?~)  and  vr  —  vl  — +  6’(7.7fl+).  In  thisca.se,  the  solution  involves 
a  leftward  moving  shock  wave  and  again  two  phase  boundaries. 

7.  Summary  of  ail  solutions.  Non-uniqueness.  In  the  preceding  section  we 
constructed  all  solutions  to  a  Riemann  problem  and  it  is  useful  to  examine  these  solutions 
on  the  {1I,vr  —  t?i,)-plane.  Note  that  7  is  the  final  strain  in  the  bar  at  large  time  while 
vr  —  vl  is  part  of  the  given  initial  data.  For  solutions  that  involve  a  phase  change,  we  map 
the  regions  /)i,  D2  and  D3  of  the  (7,7)-plane  into  the  respective  regions  D\,  and  D'^ 
in  the  {1,vr  -  n£,)-plane  by  using  the  mapping  vr  —  vi  —  G’(7,7),  7  =  7.  The  resulting 
regions  are  shown  hatched  in  Figure  8.  The  boundary  curves  A'D',  D'C',,  A' B\  I'J'  and 
K'L'  are  the  images  of  the  respective  curves  AD,  DC,  AB,  IJ  and  KL,  and  are  define  by 
VR-VL-  G{1,T{1)),  Vp~VL  =  G'(7,F’(7)),  Vr  -  Vi  =  G{%Qil)],  VR  -  VL  =  G’(7,7/,) 
and  Vr  —  vi~  G(7,7h),  respectively. 

The  solution  corresponding  to  a  point  in  D\  involves  two  phase  boundaries  and  two 
wave  fans  (Figure  6(a)).  At  a  point  in  D2  the  solution  involves  two  phase  boundaries,  one 
wave  fan  and  one  shock  wave  (Figure  6(b)),  while  in  D\  it  involves  two  phase  boundaries 
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and  two  shock  waves  (Figure  6(c)).  It  is  useful  to  note  that  the  solution  corresponding 
to  a  point  on  the  respective  curves  A'D',  D'C  and  A' B*  involves  phase  boundaries  that 
are  stationary,  sonic  and  dissipation-free  respectively. 

Solutions  which  do  not  involve  a  phase  change  can  also  be  described  on  this  plane  by 
plotting,  for  each  type  of  solution,  the  curve  vr  -  vl  =  H{l).  The  result  is  the  curve 
P'Q'R' S’  shown  in  Figure  8. 

Each  of  the  curves  P'Q’R' S',  A’B' ,A’D' ,D’C’ J'J’  and  K' L'  can  be  shown  to  be 
monotonically  increa.sing.  The  figure  reveals  two  distinct  types  of  non-uniqueness.  For 
all  initial  data  such  that  vr  —  vi  <  Hm  there  is  a  unique  admissible  solution  that  involves 
no  phase  change.  On  the  other  hand,  when  the  initial  data  is  such  that  vr  —  ujr,  >  Go  = 
G'(7o1)  703)1  there  is  a  one-pareimeter  family  of  solutions  all  of  which  involve  a  phase 
change;  the  lack  of  uniqueness  here  arises  because  of  the  undetermined  velocities  of  the 
phase  boundaries.  For  initial  data  on  the  intermediate  interval  Hm  >  vr  —  >  Gq 

both  of  these  types  of  solutions  are  available,  in  one  of  which  the  bar  ultimately  changes 
phase,  while  in  the  other  it  does  not. 

8.  Kinetics  and  nucleation.  Uniqueness.  The  propagation  of  a  phase  bound¬ 
ary  is  controlled  by  a  kinetic  relation  which,  in  the  simplest  models,  is  a  constitutively 
prescribed  relation  between  the  driving  force  /  acting  on  the  phase  boundary  and  its 
propagation  speed  i: 

fit)  =  mt))-  (69) 

The  continuum  theory  does  not  provide  an  explicit  expression  for  0;  this  is  obtained 
through  suitable  micro-mechanical  modeling.  It  is  sufficient  for  our  purposes  to  merely 
assume  that  </>  is  a  continuous  function  on  (— Ccxj,Coo)  that  increases  with  s  and  that 
the  graph  of  /  =  <^(i)  is  a  curve  that  lies  in  the  region  Fj  U  Fg  described  previously  of 
the  (s,/)-plane;  see  Figure  3.  The  entropy  inequality  (31)  implies  that  <f>(s)s  >  0  and, 
assuming  0  to  be  continuous,  that  <^(0)  =  0. 

The  kinetic  relation  controls  the  propagation  of  an  existing  phase  boundary.  A  nucle¬ 
ation  criterion  is  needed  in  order  to  signal  the  onset  of  a  phase  transformation  when  the 
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bar  involves  only  a  single  phase.  We  assume  that  the  transformation  from  the  low-strain 
phase  to  the  high-strain  phase  is  nucleated  whenever  there  exits  a  solution,  corresponding 
to  the  given  data,  with  the  associated  driving  force  /  at  least  as  great  as  a  given  critical 
value  fcr\  the  parameter  fcr  >  0  is  also  determined  by  the  constitutive  details  of  the 
material.  Once  the  nucleation  criterion  indicates  that  w^e  should  consider  solutions  that 
do  change  phase,  the  kinetic  relation  selects  the  particular  one  ol  these  solutions  that  is 
appropriate. 

The  nucleation  criterion  /(7,7)  =  fcr  describes  a  curve  M  in  the  7,7-plane  that  lies 
in  the  hatched  portion  Fa  in  Figure  7.  For  clarity  of  the  figure,  this  curve  shall  not  be 
displayed;  it  can  be  readily  verified  that  it  has  one  end  at  a  point  N  on  AD,  and  declines 
monotonically  as  7  increases  to  7— >  /?oo  as  7— >  oo.  The  coordinates  of  the  point  N  are 
found  by  solving  the  pair  of  equations 

nx^)  =  u,  =  (70) 

let  {7cr,7cr)  be  this  point  and  let  ^^(7ct,7^)  be  the  corresponding  value  of  Ct  at 

this  point.  The  image  of  the  curve  Af  in  the  (vh  —  vi,  7)-plane  is  a  curve  in  the 
hatched  region  shown  in  Figure  8  that  rises  monotonically  commencing  at  a  point  N'  on 
A'D' .  Again,  for  clarity  of  the  figure,  this  image  A^'  shall  not  be  displayed.  The  ordinate 
of  the  point  N'  is  Go--  Thus  according  to  the  nucleation  criterion,  a  solution  with  a  phase 
change  is  selected  whenever  vr  —  Vl  >  Gcr,  a  solution  without  phase  change  is  selected 
when  Vr  —  Vl  <  Gcr- 

For  vr  —  vl  >  Gcri  we  must  consider  the  1-parameter  family  of  solutions  that  involve 
a  phase  change.  In  order  to  select  the  relevant  solution  from  among  this  family,  we 
begin  by  applying  the  kinetic  relation  (69)  to  the  phase  boundaries  x  =  it  in  each  of  the 
solutions  in  Figure  6(a), (b)  and  (c);  kinetic  relation  at  the  other  phase  boundary  x  =  —st 
holds  automatically.  On  using  (7)  and  (14)  in  the  kinetic  relation  (69)  and  ensuring  that 
i  >  0,  leads  to 


iiil; 


p(7-7) 


7  =  A'(7)  for  7  >  7o3- 


(71) 
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This  describes  a  curve  K  in  the  7,^-plane  that  lies  in  the  region  Fa  (the  union  of  the 
hatched  portions  in  Figure  7);  it  commences  at  the  point  A  and  declines  monotonically 
with  /\(7)  —*  Rcc  as  7— >  oo.  When  mapped  from  here  to  the  7)-plane.  this 

yields  the  curve  K.'  defined  by  vn  —  i'l  =  G(7,A'(7)),  7  >  7o3<  see  Figure  8.  This 
curve  commences  at  the  point  A!  and  the  monotonicity  of  the  kinetic  response  function 
(p  ensures  that  it  rises  monotonically  without  bound. 

Thus  since  K,'  increases  monotonically  for  any  given  v/?  —  vl  >  Gct,  there  is  a  unique 
strain  7(>7cr)  which  satisfies  the  kinetic  relation. 

Thus,  in  summary,  we  have  shown  that  for  the  metastable  initial  data  (21).  (36).  the 
Riemann  problem  has  a  unique  solution  that  is  consistent  with  the  kinetic  relation  and 
the  nucleation  criterion.  For  initial  data  such  that  vr  —  vi  <  Gcr  this  solution  involves  no 
phase  change  and  has  the  form  given  by  Figures  5;  the  bar  commences  in  the  low-strain 
phase  and  remains  there  for  all  time.  For  vr  —  vl  >  Gcr  the  solution  involves  a  phase 
change  and  has  the  form  given  by  Figures  6;  in  this  case  the  bar,  which  was  initially  in 
the  low-strain  phase,  transforms  eventually  to  the  high  strain  phase.  In  either  case,  the 
eventual  strain  in  the  bar  is  7  and  can  read  off  Figure  8. 

In  closing,  we  note  that  some  solutions  to  the  present  problem  involved  phase  bound¬ 
aries  that  were  subsonic  while  others  involved  sonic  phase  boundaries  (the  latter  cor¬ 
respond  to  points  on  D'C  in  Figure  8).  In  the  present  analysis  the  kinetic  relation 
was  imposed  on  both  of  these  types  of  phase  boundaries  and  this  leads  to  uniqueness. 
If  the  kinetic  law  had  been  applied  only  to  subsonic  ph^lse  boundaries,  then  solutions 
corresponding  to  points  on  both  curves  fC'  and  D'C  would  be  allowed.  In  this  case 
the  Riemann  problem  corresponding  to  given  data  would  have  a  unique  solution  for 
Vr  -  Vl  <  Hm  (see  Figure  8)  but  for  vr  -  vi  >  Hm  the  problem  would  have  two  solu¬ 
tions,  one  satisfying  the  kinetic  relation  and  the  other  involving  a  sonic  phase  boundary. 
This  was  first  observed,  in  a  particular  setting,  by  Truskinovsky  [18];  he  takes  the  view 
that  the  kinetic  relation  should  only  be  applied  to  subsonic  pheise  boundaries  and  sug¬ 
gests  that  the  accompanying  non-uniqueness  (for  vr  —  vl  >  Hm  in  the  present  setting) 
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describes  an  iristabiiitv. 
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Figure  1.  Stress-strain  curve. 
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Figure  3.  Admissible  images  of  Fj  in  (5,/)-plane. 
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Figure  4.  General  form  of  solutions  to  the  Riemamn  Problem. 


1  =  low  strain  phase 


3  =  high  strain  phase 


t 


(a)  Two  shock  waves  (b)  One  shock  wave  and  one  wave  fan 


(c)  Two  wave  fans  (d)  One  wave  fan  and  one  shock  wave 


Figure  5.  Form  of  solutions  to  Riemann  problem  without  phase  change. 


Figures.  The  (7,  y/?-yr)-plane.  The  strain  7  is  the  final  strain 
in  the  bar  at  large  times,  vr  —  vl  is  part  of  the  initial 
data.  The  solution  forms  at  points  in  D\,D'2  and 
£?3  are,  respectively,  as  described  in  Figure  6(a),  (b) 
and  (c).  Points  cn  P'S'  describe  solutions  that  do 
not  change  phcLse. 
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